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We analyze theoretically the resonant inelastic x-ray scattering (RIXS) at the Ni K edge 
in the stripe-ordered state of La2-xSra;Ni04 at x = 1/3. In the calculation of RIXS spectra, 
the stripe-ordered ground state is described within the Hartree-Fock approximation by using 
a realistic tight-binding model for Ni3(i7 and 02px.y orbitals, and the electron correlations 
in the electronic excitation processes are taken into account within the random-phase ap- 
proximation. The calculated RIXS spectrum shows a tail toward the low-energy region when 
the momentum transfer of photons equals the stripe vector Q, being consistent with a recent 
experimental result. The origin of this anomalous momentum dependence of RIXS spectra 
is discussed microscopically. 
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1. Introduction 

Resonant inelastic x-ray scattering (RIXS) has been developed to be a powerful method 
for measuring elementary excitations in solids. ^"^^ This is largely owing to both the recent 
advances in instrumentations and the achievement of highly brilliant lights obtained from 
advanced synchrotron facilities. In RIXS processes, the incident photon energy is tuned to 
match an absorption energy of one of the constituent elements. The electronic system of 
materials is resonantly excited by absorbing incident photons, and then after a short time 
(typically, of the order of femtoseconds), photons are emitted, whose momentum and energy 
generally differ from those of the incident photons. The electronic system remains an excited 
state still after the photon is emitted. The momentum change and energy loss of photons 
necessarily equal the momentum and energy spent for exciting the electronic system, due 
to the energy and momentum conservation laws. Therefore, one can get information on the 
electronic excitations of the system by measuring systematically the momentum change and 
energy loss of photons. 

Excitation processes involved in RIXS depend on the material and incident x-ray 
wavenumber utilized in the RIXS measurements. Among them, RIXS at the transition-metal 
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K edges has attracted much interest, partly because the RIXS spectra arc naturally expected 
to reflect the strong electron correlations of d electrons in transition-metal compounds, such 
as cuprates,^"^^) manganites,^^'-'^'^) nickelatcs-'^^' etc. In the RIXS at the transition-metal 
K edges, transition-metal Is core electrons are resonantly excited to the transition-metal 
unoccupied bands. In the intermediate states, the transition-metal d electrons are excited 
to screen the created Is core hole. Here we should note that the core hole potential, i.e., 
the Coulomb interaction between the transition-metal Is and d orbitals, plays an essential 
role. In the final state, the 4p electron is annihilated together with the Is hole by emitting 
a photon, before the d electrons decay from the excited state to the ground state. Overall, 
the momentum change and energy loss of photons are transferred indirectly to the excited d 
electrons. 

The x-ray at the transition-metal K edges is situated in the hard x-ray regime, and 
its wavenumber is appropriate for sweeping the whole Brillouin zone. Taking advantage of 
this point, momentum dependence of the RIXS spectra has been indeed observed in various 
transition-metal compounds, such as copper oxides,^""*^^'^^'^^) NiO^'^-' etc. As a recent result 
from intensive experimental and theoretical researches, it has now become clear that the RIXS 
spectra at the transition-metal K edges refiect the charge correlation function of the strongly 
correlated d electrons. As well known, neutron scattering intensity is related to the 
spin correlation function. Thus, roughly speaking, RIXS in transition-metal compounds is to 
charge correlations what neutron scattering is to spin correlations. According to more recent 
researches, multi-niagnon excitations have been observed at low-energy regions in insulating 
cuprates.2^'2^' Therefore RIXS may provide a promising way of studying not only charge 
excitations but also magnetic excitations of strongly correlated electrons in transition-metal 
compounds in the future. 

Stripe ordering has been one of the central issues in the physics of strongly correlated 
electron systems. Neutron scattering measurements have revealed that the holes and spins 
exhibit spatial disproportionation with a stripe form in several transition-metal oxides, e.g., 
La2-a;Sr^Ni04,2^-^^) La2Ni04+5,-'^^'^^) La2-a:(Ba, Sr)3;Cu04 at x = 1,^^) etc. To try to elu- 
cidate relations to the high-Tc cuprate superconductivity, a lot of research works on stripe 
ordering have already been accumulated so far.^^^ Recently, Wakimoto and collaborators re- 
ported RIXS in two stripe-ordered 214 compounds La5/3Sri/3Ni04 and 1/8 doped La2-a;(Ba 
or Sr)a;Cu04 at the Cu and Ni K edges, respectively.^^) They observed low-energy (~ 1 eV) 
excitations with a momentum transfer corresponding to the charge stripe spatial period in the 
both compounds. The aim of our present work is to analyze this anomalous spectral feature, 
and discuss its microscopic origin for the case of La5/3Sri/3Ni04. 

The present article is constructed as follows. In § 2, we present a model Hamiltonian, 
theoretical description of the stripe-ordered state, and the formula for RIXS intensity. We 
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take the Hartree-Fock (HF) approximation for describing the stripe-ordered ground states, 
and the random-phase approximation (RPA) for taking account of the electron correlations in 
the intermediate states of the excitation processes. In § 3, numerical results of RIXS spectra 
are presented. There we present not only the results for stripe states (§ 3.2) but also for the 
undoped antiferromagnetic insulating state (§ 3.1). In § 4, the article is concluded with some 
discussions and remarks. 

2. Formulation 

2. 1 Model and Hartree-Fock approximation for stripe states 

The electronic properties of La2-xSrj;Ni04 are considered to be dominated by those of the 
Ni02 layers. Therefore we can use two-dimensional tight-binding model for a single NiO layer 
to reproduce the in-plane electronic structure of La2-a;Srj;Ni04. Since only Ni3(i7 orbitals are 
the most essential among the five Ni3c? orbitals, we take the following four orbitals: Ni3c?^2_2^2, 
Ni3d3;j2_r2, 02px and 02py. Hereafter, we specify each of these four orbitals by index £: ^ = 1 
for Ni3dj.2_2^2, £ = 2 for W\^d^^2_j.2, ^ = 3 for 02p^, £ = 4 for 02py. The noninteracting part 
of the Hamiltonian is given by 

i la i,j ll'a 

where rjj = ri — rj, £ and £' are orbital indices, cua- and ct^^ are respectively the annihilation 
and creation operators for the electron with spin a on orbital i at site i. 

For the one-particle energies, we take ei = £2 = = —9 eV, £3 = £4 = Sp = eV. For 
the hopping parameters we take ti,3(f ) = — ii,4(|) = t^p = —1.2 eV and ^2,3(1) = ^2,4(2) = 
—tdp/V^ for nearest-neighbor Ni-0 bonds, t3,4(^^) = — 13,4(^^) = tpp = 0.6 eV for nearest- 
neighbor 0-0 bonds, {x and y are the unit lattice vectors connecting inplane nearest-neighbor 
Ni sites.) The hopping matrix t(,^ii{r) satisfies the relation t(,^ti{r) = te'^i{—r) = SiSiit(^e'{—r), 
where the factor S£ equals -|-1 for £ = 1 and 2 (NiSd orbitals), and —1 for £ = 3 and 4 {02p 
orbitals). The Fourier transform of the above non-interacting Hamiltonian is given in the 
following form: 

Ho = ^^Yl (^)cLaCfe€V, (2) 

fe U' a 

with 

e«'(fc) = £e6w + te^eir), (3) 

r 

where fe denotes momentum in the first Brillouin zone (BZ). For the convenience in the 
following, we introduce the annihilation operators d and p [creation operators d^ and p^] by 

idiia,di2a,Pi3a,Pi4a) = {cn^ , Ci2a , Ci3a , Ci4a) [ 42a'PL'-PL) = (4la' 42a' 43^' 44a) ]• 

Namely, dua = when i is on Ni site and £ = 1 or 2, and pua = Cj^o- when i is on O site 
and £ = 3 or 4. For the momentum representation, we have dkea = Ckla for ^ = 1 and 2, and 
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Pkia = CMa for £ = 3 and 4. 

For the interacting part, we take the on-site Coulomb interaction at Ni sites: 

+ ^ XI XI XI AlA'a'^ii'^'^^l^''^ 
i e^i' a,a' 

+ y X X X 4iAea''^ii'^"^ii''r^ (4) 

where i is on Ni sites, i and i' are 1 or 2. [/ is the intra-orbital Coulomb repulsion, U' is the 
inter-orbital Coulomb repulsion, J is the Hund's rule coupling and J' is the inter-orbital pair- 
hopping term. Throughout the present study, we take i7 = 8 eV, U' = Q eV, and J = J' = 1 
eV. 

To deal with the many-body problem originating from H', we adopt the HF approximation. 
A lot of HF calculations using Hubbard-type models have been performed to study stripe 
ordering, so far.^^"^^) Here we introduce the mean fields 

{4eadie<7) = ^{nu + rnua), (5) 

where the nj£ and rrin are the electron number and the spin magnetic moment in orbital i at 
Ni site i, respectively. Wc assume that the stripe ordering characterized by stripe vectors Q's 
is realized. In this case, the mean fields are expressed in the form: 



me = ^e'^^'^rnQt, (7) 
Q 

where i is on Ni sites and £ is for Ni3d orbitals, and Q's characterize the spatial periodicity 
of the stripe. In the present study, wc restrict ourselves to diagonal stripe ordering expected 
from experimental results for La5/3Srx/3Ni04 and the above summation in Q is performed 
for Q = n • (n = 0, 1, 2), where = (|7r, |7r). rig^ = n*_Q^ and niQ^ = m*_Q^ hold, since 
and mn are real quantities. As a result of the HF approximation, the total mean-field 
Hamiltonian Hmf = Hq + is given in the form: 

HmF = X X ^U'i^)^iea^ke'a + X X X ^Qi<yAia^k-Qi,j + Eq, (8) 
k ll'a k Q la 



with 



^Qta = ^{nQi-mQ^a) + U' ^ nqe 
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Q e Q ti^i' 

+"1" I^ Z]("Q^"Q^' + ^Qe^he)- (10) 

Q ej^e' 

The first BZ is folded, since the spatial periodicity of the charge disproportionation is integer 
times of the original lattice periodicity. We express momentum k in the original BZ by using 
the stripe vector Q and reduced momentum fco as A; = feg + Q, where fco is restricted in 
the folded BZ. In the present study, we use the reduced BZ depicted in Fig. 1, to solve the 
HF equation for the stripe-ordered states with Qg = (|7r, |7r). Thus the total mean-field 
Hamiltonian is expressed in the form: 

Hmf = (^0 + '3)4oH-Q^aCfeo+Q^v 

fco Q U'a 

'^YYY ''^Q-Q't'jdio+Qeadko+Q'ea + ^0- (11) 

fco Q,Q' ^o- 

We should note that the summation in momentum fco is restricted only over the reduced 
BZ. By diagonalizing this mean-field Hamiltonian, we have the energy dispersions Ea^aiko) 
and unitary matrix UQ/^a,ijiko) for diagonalization, where a is the index for the diagonalized 
bands. The self-consistency condition for mean fields is given by 

= ^Yjl2J2^Q'Aa,ai'^o)^Q'+QAaA'^o)naAko), (12) 

feo Q' 
feo Q' 

where na,a{kci) is the electron occupation number on band a at momentum fco with spin a: 
'na,a{ko) = f{Ea,a{ko)) {f{E) is the Fermi distribution function). 

Here we take account of the electron-lattice interaction, since it is natural to consider that 
lattice distortions occur cooperatively with the charge disproportionation accompanying the 
stripe ordering, as discussed in refs. 41, 42 and 47. We denote the atom displacement at site i 
by u{r-i). If we assume t^^£'{r) oc |r|~"«' {aw = a^'g), the change of the hopping parameters 
due to the atom displacement is evaluated approximately by St(^(i{rij) = — Q!«'^£/'(^ij)"i7~jr-! 
where Uij = u{ri) — u{rj). Hereafter, for Uij and rij, we use units normalized by the lattice 
constant. The modification to the Hamiltonian Hq due to the lattice distortions is given by 



Throughout the present study, we take a^' = "dp = 3.5 for {£,£') = (1 or 2, 3 or 4), (3 or 4, 1 
or 2) (Ni-0 bonds), and aee = Upp = 2 for {£,£') = (3 or 4,3 or 4) (0-0 bonds) .^«) We use 
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Fig. 1. A typical reduced BZ for the diagonal stripe state is enclosed by the thick solid line. The 
wave vectors at some symmetry points are also indicated; A(-^7r, |7r), B(|7r, ^tt), C(|7r, -^tt), 
D(|^,|7r),E(f^, It:). 



the following form for the lattice elastic energy, 

EL[u^j] = ^ A'^,,A,(|r-.,|)|ti.,f , (15) 

where Ai and Aj denote the atoms (Ni or O) at sites i and j, respectively, and KAi,Aj{\fij\) is 
the elastic constant of the bond between atoms Ai and Aj. The summation with means 
that with respect to all pairs of sites. Minimizing the total energy 

^tot = (^0 + H' + Fi.d. + Ei^) (16) 

with respect to the lattice distortion Uij and using the Hellmann-Feynman theorem, we have 
the following equation for determining Uij: 

where the summation in i (i') is restricted to the orbitals on atom Ai (Aj) at site i (j). As can 
be seen easily, this condition eq. (17) is not effective for pairs of sites i and j for which t£^ii{rij) 
equals zero for any choice of i and i' . In the present study, since we take only hoppings for 
nearest-neighbor Ni3(i-02p and 02p-02p bonds, only Ujj's for these bonds are determined by 
eq. (17), while Uij^s for the other bonds are not. 

Here we assume that the lattice distortions possess the same spatial periodicity as the 
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stripe order, and then can expand Uij in the form 

UiJ = E «QW^'(ri,)e'(«-^'-^'"-^\ (18) 
Q,Q' 

where i and i' denote the atomic orbitals at sites i and j, respectively. -ffLd. is expressed in 
the momentum representation as 

^l-d- = E E E W(^0; Q, Q'yi^+Qea^^ko+Q'i'a, (19) 
fco Q,Q' U'a 

with 

w(fco; Q, Q') = -Oiw E *«.^'(^)~|7j2~ ■ '"Q-Q"AQ'-Q"^'(^)]^~'*^"'''- (2^) 

r ' ' Q" 

The energy minimization condition eq. (17) is given in the momentum representation by 

UQe,Q'e'irij) = J7 TlTTT E (^rnm'tm,m'{rij)- 

J^Ai,AA\rij\) I 



X^EE[^""'°"''^-Q,-,a,<^(^0)C/-Q',„.',„,.(fc0) 
feo ao- 

+e'''°"^^t^a',m',a,.(^o)C/Q,m,a,a(fco)]n„,.(fco), (21) 

where Ai and denote atoms on which orbitals £ and i' are, respectively, and the summation 
in m and m' is restricted to the orbitals on the atoms and Aj, respectively. UQ^£^a,a{ko) is 
the diagonalization matrix for i?MF + -f^i.d. • As easily seen, UQ£ Q'(i{rij) satisfies the following 
equalities, 

'^Qe,Q'e'{fij) = ■"■-Q^,-Q'£'(^ij) 

= -'^-Q'i'-Qei-rij)- (22) 

If we do not assume any lattice distortions, self-consistent solutions are determined by 
diagonalizing Hmf in eq. (11) and using the self-consistent conditions eqs. (12) and (13). If we 
assume lattice distortions, self-consistent solutions are determined by diagonalizing Hmf+Hi,^, 
from eqs. (11) and (19) and using the self-consistency conditions eqs. (12), (13) and (21). 

For the undoped case, we expect the checkerboard-type antiferromagnetic ground state. 
Such a ground state is obtained within the same formulation by taking Qg = (tt, tt) and 
Q = n - Qg (n = 0, 1). In this case, the reduced BZ is a square whose corners are (ibTr, 0) and 
(0, ibTr), in the wave vector space. 

We should note that the identity 

Uij + Ujk + Uki = (23) 

should always hold for any choice of lattice sites i, j and k, as easily shown from the definition 
Uij = u{ri) — u{rj). However, this identity and the energy minimization condition eq. (17) are 
not necessarily consistent with each other in general. Rigorously speaking, we should minimize 
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the total energy eq. (16) under the eondition eq. (23), but this is a rather complex and difficult 
task. In § 3.2, we present a possible resolution for this difficulty. 

2.2 RIXS formula 

To calculate the RIXS intensity, we use the useful analytic formula previously presented 
by Nomura and Igarashi.^^'^^'^^^ This analytic formula has been applied to insulating copper 
oxides,^^'^^'^^'^°) NiO,^°) LaMnOa,^-^) La2Ni04,^^) using realistic electronic structures, and 
explained the shape and momentum-transfer dependence of RIXS charge excitation spectra 
semiquantitatively. Here we outline the derivation of their formula. 

The Hamiltonian for the interaction between x-ray and electrons is given by 

= ^^^(g,e)aqe + H.c., (24) 

q,e 

Hx{q,e) = ^w{q,e)pl^^^Ska- (25) 
k,(7 

where Skcr and p\^^ are the annihilation operators for the transition-metal Is and Ap electrons 
with momentum k and spin a, age is the annihilation operator for x-ray photons with wave 
vector q and polarization vector e. The matrix element w{q, e) is given by 

p Iitt 

wtq e) = --J-^ (4p|e'«-^e • p\ Is) , (26) 
m U |g| 

in units of c = ^ = 1 (c: light velocity, h: Planck constant divided by 27r) , where e and m are 
charge and mass of the electron, respectively. 

Following Nozieres and Abrahams,^^) we calculate the inelastic scattering intensity. We 
denote the initial ground state of the electronic system without any photon by |0) in the infinite 
past time t = —oo. We assume that this state |0) absorbs an incident photon (momentum q^, 
energy Ui and polarization Cj), and emits a photon (momentum qj, energy and polarization 
e/), before a time to, and the electronic system is in an excited state |V'(*o)) at to. The 
amplitude for the state \ip{to)) is calculated within the second order perturbation theory in 



/to ru 
du dtK{to,u)HxK{u,t)HxK{t,-oo)al^^\0), (27) 
-oo J —oo 

where K{u, t) is the time evolution operator in the case of Hx = 0, and t and u are the times of 
the photon absorption and emission, respectively. The total number of x-ray photons generated 
before to with wave vector qj^ and polarization e/ is Nf{to) = ('0(io)|"qj-e/aq^e/|'0(*o))- 
Inelastic x-ray scattering spectra are regarded as the number of photons generated in a unit 
time. Thus, we obtain the scattering intensity by deriving Nf{to) with respect to to and 
contracting photon annihilation and creation operators: 



/oo f'u' f'O 

du' / dt' / dte'^*^^'-^^e-''^f''' S{t,0;u' ,t'), 
-oo J —oo J —oo 



(28) 
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with 

S{t,u-u',t') = {0\Hl{t'-qi,ei)H^{u'-qf,ef)Hl{u-qf,ef)H^^^^^ (29) 

where Hx{t;q,e) is the interaction representation of Hx{q,e). The function S{t,u;u' ,t') can 
be calculated by the Keldysh formalism. 

For the case of RIXS at the transition-metal K edges, we introduce the following ap- 
proximations: (i) We take completely flat dispersion for the Is band, since the Is electrons 
are strongly localized in the inner K shell, (ii) We use a free-electron model for the 4p elec- 
trons, since the transition-metal 4p orbitals are expected to much extend in space. In the 
present study, we take simply a cosine-shaped band for the 4p electrons. This is justified by 
the fact that the excited 4p electron plays only a role of "spectator" as far as we discuss 
the momentum dependence of RIXS spectra. Of course, detailed 4p electron energy disper- 
sions are necessary for more quantitatively precise discussions on K edge absorption spectra 
and dependences on the incident photon energy and polarization, (iii) Since transition-metal 
d electrons also possess a localized nature, the Is core-hole potential (whose absolute value 
equals the Coulomb integral Vis-d between the Is and d electrons) is expected to be rather 
strong. Nevertheless, we use the Born approximation for Vig-d- The Born approximation was 
partly justified by taking account of multiple-scattering processes for the case of La2Cu04.22) 
On the basis of these assumptions, we can perform the integrals with respect to time variables 
in eq. (28), and have 



W{qiei;qfef) = 4 

with 



Vis-d ^ w{qi,ei)w{qf,ef)* 



E 



N ^ -i{ui\k)-i{uf\k) 



2 



7(w;fe) = a; + £is + iris-£4p(fe), (31) 

where £is. Tig and £4p(fc) are the Is energy level, the decay rate of the Is core hole, and the 4p 
band energy, respectively. Y4{q,Lj) is the Fourier transform of the dynamical charge-density 
correlation function for d electrons: 

Yd{q,t) = {pdq{t)Pd-qm, (32) 

where Pdq{t) is the Heisenberg representation of the following electron density operator for d 
electrons, 

k,cr 

We should note that the Is core hole plays a role of a localized testing charge perturbing the 
d electrons through the Is-d Coulomb interaction. As easily seen from the above discussions, 
applying the Born approximation to Fis-rf is equivalent to taking account only of the linear 
response of d electrons to the perturbation due to the Is core hole charge. The fluctuation 
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dissipation theorem relates Yd{q,uj) to the charge susceptibihty of d electrons XdiQ,^)'- 

Ydiq,Lj) = 2(1 - e-'^/^)-ilmxd(q,a; + iO). (34) 

In the present study, the incident photon energy Ui is set near the K absorption energy, 
i.e., cji £4p{kedge) ~ ^Ti^. Then 7(0;^; k) becomes small and the RIXS intensity W{qiei; QfBf) 
is resonantly enhanced. The resolution of the RIXS spectra as a function of photon energy 
loss oj = uji — ujf is determined by that of Yd{q,uj), i.e., the decay rate of charge excitations 
of d electrons, rather than by the core hole decay rate Fi^. 

To proceed further, we have to evaluate explicitly Yd{q,uj). Specifically we take the HF 
approximation to describe the stripe-ordered ground state of the d electrons, and the RPA to 
take account of the electron correlations in the intermediate state. RPA means that we neglect 
any couplings among various modes corresponding to different wave vectors qr's, in other words, 
we assume various excitation modes with different qr's are separately renormalized by electron 
correlations. Consequently we have the final expression for stripe-ordered states: 

W{qiei,qfef) = 8n ^ 



E 



^ 7(a;j;fc)7(a;/;fc) 



h^^^i^ + EaA^o) - -Ea/,(,({fco + Q}))na,a(fco)[l - «a',a({feo + Q})] 



X 

N 

fco aa' <j 



W Q,Q' 



(35) 



where Qk^+q is the stripe vector by which the wave vector fco + Q is pulled back into the 
reduced BZ, and then the wave vector fco + Q is reduced to the wave vector {fco + q} in the 
reduced BZ, i.e., {fco -|- q} = fco + q — Qk^+q- ^e'eiQ) is the vertex function including electron 
correlation effects, which we calculate within RPA with respect to the Coulomb interaction H' 
of eq. (4). On the other hand, if we turn off the RPA corrections by setting AJ,^{Q) = 5q5^i(^, 
then we can extract simple band-to-band charge excitations. 

Numerical diagonalization technique using finite cluster models is another promising 
method for calculating RIXS intensity, in the sense that it can take full account of 
electron correlations without any approximations, but may have severe difficulty in analyzing 
detailed momentum dependences for such stripe-ordered systems as La5/3Sriy3Ni04, since the 
unit cell in stripe states contains a relatively large number of atoms. 

3. Numerical Results 

3.1 Case of undoped antiferromagnetic insulating state: x = 0, Q = (tt, tt) 

Before proceeding to the cases for stripe ordered states, we present the results for the 
undoped antiferromagnetic state. Based on the formulation in § 2, the antiferromagnetic 
ground state is obtained: the electron occupation number and total staggered spin moment 
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are rii = 2.24 and mj = 1.74 (in units of /^b) at each Ni site i, respectively, where we have 
not assumed lattice distortions. The calculated RIXS spectra for three momentum transfers 
q = (0, 0), (tt, 0) and (tt, tt) are presented in Fig. 2. We find three spectral features around the 
energy loss w = 4, 6 and 7.5 cV. Overall, the spectral peak positions in the calculation seem 
consistent with the experimental results semiquantitatively. On the other hand, the total in- 
tegrated intensity at q = (0, 0) seems clearly larger than that at q = (vr, 0) in the experiment, 
while it does not seem in the theory. This inconsistency may be resolved by taking account 
of multiple scatterings beyond the Born approximation for the core-hole potential, since the 
low-energy intensity at q = (0, 0) seems more enhanced due to the multiple scatterings than 
at q = (vr, 0), according to the study for La2Cu04.^^) In addition, for more complete quan- 
titative consistency in spectral shape, we should use not such simple tight-binding electronic 
structures but more precise electronic band structures for all of the Ni3d, 02p and Ni4j) states. 
Comparing the peak positions between the cases with and without RPA, the 4 eV peak is 
shifted to the high energy region due to electron correlations at q = (vr, tt), while it is not at 
q = (0,0) and (7r,0). 

Here, we present more detailed momentum dependence of the RIXS spectra than presented 
in the previous work by Takahashi et al.^^^ The detailed momentum dependence of RIXS 
spectra along the symmetry lines in the first BZ is shown by the intensity plot in Fig. 3. Our 
calculation suggests the possibility that the 4 eV peak exhibits stronger dispersion along the 
line (0, 0)-(7r, vr) rather than along the line (0, 0)-(7r, 0). This contrasts strongly with the case 
of La2Cu04, in which the 2 cV peak shows strong dispersion along (0, 0)-(7r, 0).^^ 

To see the microscopic origin of the RIXS weights, we present the results of density of 
states (DOS) in Fig. 4. The insulating gap energy, which corresponds to the charge transfer 
energy between the Ni3d and 02j; bands, is about 4 eV, being consistent with photoemission 
experiments. ^'^^ The 4 eV RIXS spectral weight corresponds to the charge transfer excitations 
from the Ni3d-02p anti-bonding band to the NiSd upper Hubbard band. The 6 eV peak cor- 
responds to the charge transfer excitation from the Ni-0 bonding band to the upper Hubbard 
band. 

3.2 Case of doped stripe- ordered states: x = ^, Qg = (|7r, |7r) 

As mentioned in the last paragraph of § 2.1, it may be impossible to determine all w^'s 
by solving eq. (17) consistently with eq. (23). For example, let us focus on a Ni atom and its 
two neighboring O atoms whose two Ni-0 bonds cross perpendicularly each other. If those 
two Ni-0 bonds shrink (the 0-Ni-O angle remains perpendicular) , then the 0-0 bond length 
determined from eq. (23) necessarily shrinks. However, this shrinkage of the 0-0 bond is 
not necessarily consistent with Uij determined for that 0-0 bond by using the condition 
eq. (17). For this difficulty, we may take the following way: firstly we determine uij for each 
nearest-neighbor Ni-0 bond by using the condition eq. (17) (or equivalently eq. (21)), and 
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Fig. 2. (Color online) Calculated results of RIXS spectra for three momentum transfers q = 
(0, 0), (tt, 0), (tt, tt), and comparison with experimental results at q = (0, 0), (tt, 0). The empty 
circles represent experimental data read from Ref. 15. The thick red and thin blue lines represent 
the results calculated with and without RPA, respectively. 



then we determine Uij for each nearest-neighbor 0-0 bond by using eq. (23). Fortunately, 
this treatment seems valid, as inferred from the following discussions. As far as we studied, if 
we determine all Ujj's (including those for 0-0 bonds) from eq. (17) (or equivalently eq. (21)) 
without using the condition eq. (23), then the electronic structure and ■Ujj's only negligibly 
depend on i^o,0) while they depend strongly on Kni.o- This means that the lattice distortions 
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Fig. 3. (Color online) Intensity plot of the RIXS spectra calculated for the undoped antifcrromagnetic 
state. The horizontal axis represents the momentum transfer along the symmetry lines of the square 
BZ, and the vertical axis represents the energy loss of photons. 
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Fig. 4. (Color online) Total and NiSd partial density of states are depicted. The total DOS, and the 
partial DOS for ^\id,j.i^yi and Ni3(i322_r2 orbitals are represented by thin black, thick red solid 
and thick blue broken lines, respectively. 



are strongly dominated by changes of Ni-0 bond length and the contributions from 0-0 bonds 
are almost energetically negligible. Therefore, to determine Uij for each nearest-neighbor 0-0 
bond, we may neglect the condition eq. (17) for 0-0 bonds, and instead of that, we should 
use eq. (23). The condition eq. (23) for 0-0 bonds reduces in the momentum representation 
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where £ is Ni d orbital (i.e., ^ = 1 or 2), and the double signs in front of | (or |) correspond 
in each equation. 

We consider the following four self-consistent solutions (I)- (IV). (I): not allowing lattice 
distortions, we diagonalize Hmf from eq. (11) and use the self-consistency conditions eqs. (12) 
and (13). (II): allowing lattice distortions, we diagonalize i?MF + -f^Ld. from eqs. (11) and (19), 
and use the self-consistency conditions eqs. (12), (13) and (21), where we use eq. (21) only for 
Ni-0 bonds and do not use eqs. (36) and (37) for 0-0 bonds. For (II), we take KNi,o = 30 
eV. (Ill) and (IV): allowing lattice distortions, we use eq. (21) only for Ni-0 bonds, and use 
eqs. (36) and (37) for nearest-neighbor 0-0 bonds. For (III) and (IV), we take -fCNi,o = 30 eV 
and KNi,o = 60 eV, respectively. These solutions are summarized in Fig. 5. Concerning lattice 
distortions in (II), nearest-neighbor Ni-0 bonds shrink particularly around Ni sites with excess 
hole density (Q sites in Fig. 5), being qualitatively consistent with previous studies.^^'^^^ 

In Fig. 6, we present the calculated RIXS spectra in the low-energy region at various q 
points on the symmetry lines, for the stripe solution (III). A remarkable feature is that the 
low-energy edge of the spectra shows a tail toward the low-energy region at q = Q = (|7r, |7r). 
This is consistent with experimental results by Wakimoto ct al.'^'') It should be noted that this 
low-energy tail appears already in the spectrum calculated without RPA, as shown in Fig. 6. 

To see more detailed momentum dependence of the RIXS spectra on momentum transfer 
q, we show intensity plots of the calculated RIXS spectra along the diagonal lines of the 
original square BZ for the stripe solutions (I)-(IV) in Fig. 7. Comparing Fig. 7 with Fig. 3 of 
Ref. 37, we consider that the RIXS intensity for the stripe state (III) is the most consistent 
with the experimental result. Comparing the spectra for the stripe states (II) and (III), we can 
see that shrinkage of nearest-neighbor 0-0 bonds, which is caused by shrinkage of nearest- 
neighbor Ni-0 bonds, enhances the low-energy dispersive behavior around q Q. Comparing 
the spectra for the stripe states (III) and (IV), we can sec that the low-energy gap at q = Q 
becomes larger for smaller i^Ni,0) i-e., for softer Ni-0 bonds. 

So far, we have considered only a single stripe-ordered domain characterized by the stripe 
vector Qs = (|vr, |7r). To compare with the experimental results in more detail, we should 
bear in mind that, in addition to the above kind of domains, actual samples will contain also 
domains which correspond to Qg = (|7r, |7r) (sec Fig. 8(a)), as one can expect easily from 
tetragonal lattice symmetry. Here we assume that those two kinds of domains take evenly 
the same volume fraction in the sample, and the domain walls separating the domains only 
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no na 
rriQ m/\ mu 


(I): no lattice distortions 


2.16 2.20 2.16 
1.59 -1.77 1.59 


(II): i^Ni,o = 30 eV, 

not using eqs. (36) and (37) 


2.26 2.37 2.34 
0.98 -1.45 1.40 


(III): KNi,o = 30 eV, 

using eqs. (36) and (37) 


2.30 2.39 2.35 
1.06 -1.42 1.31 


(IV): K^i,o = 60 eV, 

using eqs. (36) and (37) 


2.22 2.29 2.27 
1.24 -1.64 1.50 



Fig. 5. (Color online) (a) Schematic figure of diagonal stripe states with — (|7r, |7r). The same 
symbols represent equivalent Ni sites, (b) Typical stripe solutions within the HF calculation are 
presented, n^'s and m^'s (i — 0: '^i D) are the electron occupation number and total spin moment 
(in units of /ie), respectively, at Ni sites marked with the same symbols in (a). 



negligibly affect the RIXS spectra. Under this simple assumption, the total RIXS spectrum is 
given by the average of the contributions from the two kinds of domains. The averaged RIXS 
spectra for the stripe solution (III) are presented in Figs. 8(b) and (c) (compare Fig. 8(c) with 
Fig. 3 of Ref. 37). 

According to the theoretical spectra (the upper two panels of Fig. 6 for q = (0, 0) and 
(vr, vr), and the two panels of Fig. 8(b) for q = (vr/S, vr/3), q = (27r/3, 27r/3), q = (47r/3, 27r/3) 
and (57r/3, 7r/3)), we see that the intensity at g = (0,0) and (vr, vr) seems relatively weak 
compared with that at the other momentum transfers. One might consider that the weakness 
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Fig. 6. (Color online) Calculated RIXS spectra in the low-energy region at various momentum trans- 
fers for the stripe solution (III). In each panel, the thick red and thin blue lines represent the 
results calculated with and without RPA, respectively. 



of the intensity at q = (0, 0) is inconsistent with the strength of the experimental intensity 
at the point 'A' in Fig. 2 of Ref. 37. However, we should note that, in general, the intensity 
around the F point (i.e. q = (0,0)) can be affected easily by several realistic factors which 
are not included in the present model. One of them is the effect of disorders in the sample. In 
doped systems, actual samples inevitably contain some disorders or disordered domain walls 
situated randomly, while doped carriers are assumed to form ideally periodic configuration 
in the theoretical model. In principle, such disordered fractions possessing no characteristic 
spatial periodicity contribute to the intensity at the F point. Another possible factor is the 
long-range component of Coulomb interaction. In doped systems, the long-range Coulomb 
interaction can change the form of the RIXS spectra around the F point. In addition, the 
multiple scattering due to the core-hole potential may enhance the spectral intensity around 
q = (0,0), as mentioned in § 3.1. To take account of these effects thoroughly remains an 
interesting but still difficult future work. 

To study the origin of the above-mentioned anomalous momentum dependence of RIXS 
spectra around q ~ Q, we present the electron energy dispersions for the stripe solutions in 
Fig. 9. In Fig. 9, only the majority (i.e., up) spin bands are depicted for the cases (I), (III) and 
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Fig. 7. (Color online) Panels (a), (b), (c) and (d) show the intensity plots of calculated RIXS spectra 

for the four stripe solutions (I), (II), (III) and (IV), respectively 
(see the text and Fig. 5 for the four stripe solutions). In each panel, the horizontal and vertical 
axes represent momentum transfer q and energy loss uj, respectively, and q sweeps along the 
diagonals of the square BZ. Q = (|7r, |7r) is the stripe vector. 



(IV). The calculated bands suggest the existence of gap minima around the symmetry points 
A, B and C (and the equivalent points), as seen in Fig. 9. As easily expected from this, the 
low-energy tail around q ~ Q in RIXS spectra originates from the charge excitations between 
those symmetry points. Such gap minima may be observable by angle-resolved photoemission 
spectroscopy (ARPES). 

4. Discussions and Conclusions 

We have presented the RIXS spectra calculated for undoped insulating antiferromagnet 
La2Ni04 and stripe-ordered La5/3Sri/3Ni04, where we have described the ground state by 
the HF approximation and have taken account of electron correlations within RPA. For the 
undoped case, we have explained the spectral peak positions semiquantitatively, and have pre- 
sented a more detailed plot than in the previous work. In the experiment by Collart et al.,^^^ 
the momentum-transfer dependence along the line (0, 0)-(7r, 0) was studied in detail. However, 
our present calculation suggests that the low-energy peak around 4 eV shows stronger dis- 
persion along (0, 0)-(7r, vr) rather than along (0, 0)-(7r, 0). This contrasts with the case of the 
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Fig. 8. (Color online) (a) Schematic figure of two kinds of diagonal-stripe domains, (b) Averaged 
RIXS intensity from the two stripe domains at q = (vr/S, tt/S), (27r/3, 27r/3), (47r/3, 27r/3) and 
(57r/3, 7r/3). The thick red and thin blue lines represent the results calculated with and without 
RPA, respectively. The spectra at q = (0, 0) and (tt, tt) are unchanged by taking the average, and 
are presented still by the upper two panels of Fig. 6. (c) Intensity plot of the averaged RIXS spectra. 
The horizontal and vertical axes represent momentum transfer q and energy loss w, respectively, 
and q sweeps along the diagonals of the square BZ. 

insulating cuprate La2Cu04. Experimental verification of this suggestion is awaited. 

For doped stripe-ordered La5/3Sri/3Ni04, we have explained the anomalous momentum- 
transfer dependence of spectra observed experimentally, i.e., the calculated RIXS spectra show 
a tail toward the low-energy region when the momentum transfer of photons equals stripe vec- 
tor Q = (|7r, Ivr), being consistent with the recent experimental result by Wakimoto et al.'^'^^ 
The reason for the low-energy tail at Q is that the gap minima of electron energy disper- 
sion in the stripe-ordered states exist around the symmetry points A(^7r, |7r), B(|7r, ^vr), 
C(|7r, ^vr) and the equivalent points in momentum space. This feature of the energy bands 
in stripe-ordered La5/3Sri/3Ni04 may be verified by ARPES experiments. 

We searched for other possible self-consistent solutions than presented in Fig. 5 (b), by 
choosing various initial values of nn, mn and Uij in numerical iterations. However, we found 
only the equivalent solutions, which are obtained by rearranging {nQ,mQ), {n/\,m/\) and 
{no, "ma) or by reversing the signs of spin moments. The stripe states with lattice distortions 
((II), (III) and (IV) in Fig. 5 (b)) are not symmetric under spatial inversion. This does 
not mean that the electron energy and RIXS spectra are not symmetric under inversion in 
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Fig. 9. (Color online) Energy dispersions along the symmetry lines of the reduced BZ (Fig. 1) in the 
four stripe states (I)- (IV) (see the text and Fig. 5). For (I), (III) and (IV), only up (majority) 
spin bands are depicted, while both of the up and down spin bands are depicted for (II). Each 
symbol. A, • • • , E, on the horizontal axis denotes the k point with the corresponding symbol in 
Fig. 1. The Fermi level is set to zero. 



momentum space. In fact, the obtained electron energy dispersions are symmetric (see the 
dispersions along the D-E line in Fig. 9). We checked numerically that the RIXS spectra are 
also symmetric under inversion in wave vector space, as a result from the symmetry of the 
electron energy dispersion. 

For the case of stripe ordered states, we have presented theoretical results of the RIXS 
spectra only for the relatively low-energy region (0-2. 5eV, as in Figs. 6, 7 and 8). We should 
note that the RIXS spectra calculated for doped stripe states within the present formulation 
are reliable only in such low-energy region, while those for the undoped antiferromagnetic state 
are reliable up to the high-energy region. For the undoped case, the HF approximation presents 
large magnetic moments {mn = 0.78, = 0.96), and consequently the two magnetically 
split bands with an enough large energy gap mimic well the actual Hubbard bands, which are 
separated by the Mott-Hubbard gap of the order of U. In this case, the electronic structure 
in the actual antiferromagnetic ground state of the Mott insulating phase is well described 
over a wide energy range up to the order of U by the HF approximation. On the other 
hand, in the doped stripe states, the magnetic moments are reduced {mn = 0.22-0.54, mi2 = 
0.85-0.88, for the stripe solution (HI)), compared with the undoped case. In this case, the 
electronic structure only inside the reduced magnetic gap can be still described by the HF 
approximation, although the high-energy Hubbard satellite bands are no longer reproduced. 
Very roughly speaking, the magnetic gap for orbital £ within the HF calculation is evaluated 
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to be Aj£ ~ rriii'. According to this rough cvahiation, Aj^ ~ 2.6-7.5 eV for the 

stripe solution (III), while An ~ 7.2-8.5 eV for the undoped antiferromagnetic ground state. 
Thus, for the doped stripe-ordered cases, the calculated RIXS spectra are reliable only in the 
low-energy region up to about 2.6 eV. To study RIXS spectra in such doped Mott insulators 
over a wider energy range, we require more advanced methods such as the dynamical mean- 
field theory (DMFT).^^) 

As we have discussed, the bond shrinkage between nearest-neighbor 0-0 sites seems im- 
portant to obtain that low-energy anomalous dispersion quantitatively consistent with the 
experimental data. This is related to the fact that the doped holes occupy mainly the 02p 
orbitals. The RIXS spectra reflect the momentum distribution of the partial component of 
the Ni3(i states mixed with the 02p states. Therefore, change of the 0-0 bonding length (or, 
in other words, change of the 02p band width) affects strongly the low-energy RIXS spectra, 
although the contribution of 0-0 bonding is energetically negligible as mentioned in § 3.2. 

In conclusion, we would like to point out that the anomalous low-energy RIXS weight 
at the stripe vector Q does not indicate some kinds of collective charge excitation modes, 
since it can be explained at least qualitatively within simple band-to-band transitions. This 
physical picture seems consistent with the robustness of the charge ordering observed under 
high electric fields. ^^'^'^^ 
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